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Ph! Abstract 

Oh: 

^ . We study QED corrections to chiral symmetry breaking in the Nambu-Jona- 

Lasinio (NJL) model with two flavors of quarks. In this model, the isospin symmetry 
is broken by the differences between the current quark masses and the electromagnetic 
charges of the up and down quarks. To leading order in the expansion, we calculate 
the effective potential of the model with one-loop QED corrections at finite tempera- 
ture. Evaluating the effective potential, we study the influence of the isospin symmetry 
breaking on the orientation of chiral symmetry breaking. The current quark mass plays 
an essential role in maintaining the orientation of the chiral symmetry breaking. If the 
average of the up and down quark masses is small enough, we find a phase in which 
the pion field has non-vanishing expectation value and dynamical CP violation takes 
place. 
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§1. Introduction 



Quantum chromodynamics (QCD) is the first principle in describing the physics of quarks 
and gluons. The QCD Lagrangian has the global flavor symmetry SUi,{Nf)® SUK{Nf) for 
the A^f flavors of massless quarks. In nature, quarks have different current masses for all 
flavors, and the global flavor symmetry is an approximate symmetry for light quarks. Owing 
to asymptotic freedom in the QCD interaction, the coupling constant grows at low energy 
scales, and the hadronic phase is realized. In the hadron phase, the approximate chiral 
symmetry is dynamically broken. 

The broken chiral symmetry is restored in certain extreme environments, e.g. high tem- 
perature and high density. This restoration of chiral symmetry is a non-perturbative phe- 
nomenon in QCD. The phase structure of QCD has been studied in some low energy effective 
theories^^"^^ and through numerical calculations in lattice QCD. Prom these results, it is con- 
jectured that the chiral symmetry should be restored near the critical temperature, T < 200 
MeV. Recently, an experimental study at RHIC found evidence of symmetry restoration 
from hadronic matter to a state of deconfined partonic matter. 

The effect of an electromagnetic field on chiral symmetry breaking is complicated. Low 
energy effective theories provide possible approaches for studying non-perturbative QCD phe- 
nomena in an external electromagnetic field. The NJL model is the simplest model in which 
the chiral symmetry is broken dynamically. The dynamical origin of symmetry breaking in 
the NJL model has been studied in the case that there exists an external electromagnetic 
field. The combined effect of both an external electromagnetic field and gravity has also 
been investigated in Refs. 6) and 7). Then it was found that the external electromagnetic 
fields can induce a very wide variety of phases. 

Radiative QED corrections play an important role in the physics of quarks and hadrons, 
even if we consider a system with no external electromagnetic field. The difference between 
the electromagnetic charges of up and down quarks breaks the SU{2) isospin symmetry. The 
difference between the masses of up and down quarks also breaks the isospin symmetry. 
Because of the isospin breaking, it is conjectured that a variety of phases are realized inside 
quark matter through QED corrections. If the sum of the up and down quark masses, 

+ ^id, is small enough, there is a possibility that a pion field will develop a non- vanishing 
vacuum expectation value. 

We assume that the NJL model is a phenomenological low energy effective theory of QCD 
and that the fundamental QCD Lagrangian is invariant under the CP transformation. Thus, 
the quark mass should be real and there should be no ^-term. Because of the anomaly, the 
chiral U{1) transformation, which flips the sign of the quark mass, produces the ^-term in 
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QCD. This is only a phase convention for quark fields, and such a ^-term does not induce any 
CP violating phenomena. The situation may be different after dynamical chiral symmetry 
breaking. The orientation of the vacuum is fixed, and thus the CP symmetry can also be 
broken. This is an example of the dynamical CP violation proposed by Dashen.^^ To observe 
the phase structure of the theory, it seems more convenient to choose a phase convention for 
the quark fields in which the Lagrangian has no ^^-term. There is no degree of freedom that 
allows us to set the quark mass to a positive value in this convention. For this reason, both 
positive and negative quark masses should be considered. 

In the present paper, we study the infiuence of the current quark mass and the quan- 
tum corrections of the electromagnetic fields on the chiral symmetry breaking in the NJL 
model. The order parameters of the chiral symmetry breaking are given by the vacuum 
expectation values of the a and tt" fields. These expectation values are found by deter- 
mining the minimum of the effective potential. First, we calculate the effective potential of 
the NJL model to leading order in the expansion. Evaluating it, we derive the phase 
structure of the NJL model. Next, the gauged NJL model is considered. We calculate the 
effective potential with one-loop QED corrections and evaluate the phase structure of the 
gauged NJL model. The charged and neutral pion mass difference is investigated in §4. As 
is well-known, a higher-order correction in the expansion is essential to elucidate the 
pion mass difference. We introduce it through phenomenological meson kinetic terms. The 
phase structure is evaluated in the gauged NJL model with additional meson kinetic terms. 
In §5 we study the temperature effect in the imaginary time formalism. The behavior of the 
effective potential is evaluated in both the NJL model and the gauged NJL model near the 
critical temperature, T ~ 170 MeV. In the case of small \m-a + md|, there is a parameter 
space in which the pion field acquires a non-vanishing vacuum expectation value. In §6, 
the dynamical origin of CP violation is studied within the NJL model for small |mu -|- md|. 
Finally, we give some concluding remarks. 

§2. NJL model with finite current quark meiss 

We start from the Lagrangian density of the two-flavor NJL model with flnite current 
quark mass. It is deflned by 

= i^{iYd, -M)^ + ^ [i^^f + (V^Z75r»'] , (2-1) 

where G is the coupling constant of the four-fermion interaction, A'^ represents the number 
of colors, M is a 2 X 2 matrix which contains the current quark masses for the up and 
down quarks, and r° {a — 1,2,3) are the Pauli matrices of the isospin vector. We assume 
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a diagonal form for the mass matrix, M = diag(mu, mA). In our phase convention it 
cannot be assumed that the mass matrix is positive. In four dimensions this model is not 
renormalizable, because the coupling constant G scales with the mass as (mass)~^. The 
model is defined with a regularization parameter and is regarded as a low energy effective 
theory of QCD. 

In the massless limit, M — > 0, the Lagrangian density (12 -ip is invariant under the isospin 
symmetry SU\{2) transformation, il){x) exp(i6"'T"' /2)iIj{x) , and the chiral transforma- 
tion, iIj{x) exp(z75^^5r"/2)?/'(x). In this case, the chiral charge is defined by = 
J d^x^/'7o75('r'^/2)^/'. The commutation relations between the chiral charge and fermion bi- 
linears are given for the composite scalar operator ipip as 

[zQl ^Pij{x)] = ijzj,T-ij{x) (2-2) 

and for the composite pseudo-scalar operator ■ipi'y^T^'ip as 

[tQl ^tj,T'4j{x)] = -d'^'^^ix). (2-3) 

If the composite operators tpij^r'^tlj and/or i/jip develop non- vanishing expectation values, 
the chiral symmetry is dynamically broken. 

In the real world, up and down quarks have small current masses. The mass term is not 
invariant under the chiral transformation; it breaks the chiral symmetry explicitly. In the 
case mu 7^ mj, the isospin symmetry, SUy{2), is also broken. 

For simplicity, we rewrite the Lagrangian density (12-11) using the auxiliary fields method. 
First, we introduce the auxiliary scalar field cr ^ —{G/N)ipip and pseudo-scalar fields vr" ~ 
— {G/N)ipi'y5T°'ip. Thus the Lagrangian density (12-ip is rewritten as 

= tPiiYd, -M-a- z75r'^7r")V' - ^ + (^")'] • (2-4) 
A neutral pion and charged pion fields tt^ are defined by 

= r"7r". (2-5) 
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To study the phase structure of the model, we evaluate the effective potential. The 
ground state is found by finding the minimum of the effective potential. Here, we start with 
the following generating functional of Green functions: 



= / Vi/jVi/jVaVn exp ( i / d^^xC^^ 



VaVnexp ^^^j-^ j d^x [a^ + {n^'f] 
i In det(i7''9^ -M-a- ij^r^Ti" 



[2-6) 
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Owing to translational invariance, the expectation value of the auxihary fields should be 
constant. Furthermore, the quantum corrections from the internal lines of the auxiliary 
fields disappear in the large N limit. To leading order in the expansion, the path 
integral over the auxiliary fields in Eq. fl2-6p gives only an overall factor. Thus the effective 
potential is found to be 

1 

2iNjd^x 

= ^ [0-2 + (7r«)2] + lndet(z7'^a, -M-a- t^.rV). (2-7) 

The current quark masses for up and down quarks are much smaller than the vr meson 
scale. For this reason, we evaluate the effective potential up to second order in the current 
quark mass and apply the expansion 

i\ndet{iYdf, -M-a- i-f^r^Tc") 

00 

= itr ln(i7^9^ -a - i'^^t'^'k'') + tr ^ /„ 

n=l 

= iii\Yi{iYd^ -a- t^^rV) + tr(/i + h) + ©(m^), (2-8) 



where tr represents the trace over the flavor, spinor and space-time coordinates, and is 
given by 



1 / M 



in \i7^9^ — a — i'j5T"-7r°- + ie 



(2-9) 



In four dimensions, Eq^ fl2-8p is divergent. To obtain a flnite result, we introduce the 

After some straightforward calculations, a finite expression for 



three-momentum cutoff A 



*)| 



the effective potential is obtained: 



{m^ + ma)(xg{cx'^;Af^) 



2 
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*■* We use a three-momentum cutoff to obtain a finite result, employing the same regularization method 
even at finite temperature. 



5 




-1x10-" 



-0.4 -0.2 0.2 0.4 

O(GeV) 

Fig. 1. Behavior of the effective potential y(cj, vr" = 0) for mu + nid = MeV (red curve), + 
= 4 MeV (blue curve) and + = 9 MeV (black curve) . 



where we have defined a = a + (tt") , and the functions f{s , t ) and g{s , t ) are given by 
f{s'; e) ^ {2e + s') - In "^^^^f^ , (2-11) 



g{s',e) ^ Vt^(t^ + .^)-.^ln ^+^^f + ^ . (2-12) 



The effective potential (12-101) is symmetric with respect to the isospin transformation 
up to terms hnear in the current quark mass. Because the isospin breaking contribution 
appears only at orders higher than O(m^), it gives only a small effect. For = — ^ 
the terms of 0{m) in Eq. (2.10) vanish and the terms of 0{m?) mainly contribute to the 
orientation of the chiral symmetry breaking. It induces the pion condensation. To determine 
the contribution from the isospin breaking, we evaluate the behavior of the effective potential 
(12-101) as the difference between the masses of the up and down quarks varies. In Table [U we 
list the values of a and vr at the local minimum of the effective potential in the vr" = plane, 
V^((j, vr" = 0), and the cr = plane, Vm(cr = 0,7r"), respectively. We also list the values of 
the effective potential for both points. The vacuum state corresponds to the point with the 
smallest value of the effective potential. In our model, this state is located in the a = or 
7r° = plane. 

As seen in Table HI the quark mass mainly contributes to the a direction of the effective 
potential, while the isospin breaking in the quark mass term induces the (7r^'^)o-=o,7r3=o ~ 
(7r^)o-=o,7ri.2=o difference. However, the ratio of the depths of the effective potential, 

V{a = 0, (7ri'2),=o,.3=o, vr^ = 0) - V{a = 0, tt^'^ = 0, {'k''),^o,.^.2=^) 



(2-13) 



Via = 0, (7ri.2)^^o ,^3=0, = 0) + V{a = 0, ttL^ = o, (7r3),=o,.i.2=o) 

is approximately 0.008% for 'm^^ = 2.5 MeV and ma = 6.5 MeV. 

We plot the behavior of the effective potential in the vr" = plane in Fig. [H It is seen 
that increasing the total current quark mass of the up and down quarks, + ma, causes 
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Table I. Influence of the isospin breaking on the effective potential. The columns labeled ((T)jra=o 
and (7r'*)o-=o list the local minima of Vm{(T,TT'^ = 0) and Vm{o- = 0, vr"), respectively. 



(MeV) 


(MeV) 


(GeV) 


(GeV) 


Kn(((T)^"=0,^" =0) 

(GeV*) 


K.(a = 0, (7r'^).=o) 
(GeV*) 


4.5 


4.5 


0.2829 


0.2719 


-9.107 X 10-5 


-6.610 X 10-5 


2.5 


6.5 


0.2829 


0.2718 


-9.105 X 10-5 


7ri'2 : -6.608 x 10-5, 
TT^ : -6.609 X 10^5 


-4.5 


4.5 


0.2718 


7ri^2 ^ 0.2718, 
TT^ = 0.2719 


-6.604 X 10-5 


7ri'2 : -6.604 X 10^5, 
TT^ : -6.610 X 10^5 


-2.5 


6.5 


0.2769 


0.2718 


-7.706 X 10-5 


7ri'2 : -6.603 x 10-5, 
TT^ : -6.609 X 10^5 



the effective potential to increase for a negative a and decrease for a positive a. Thus the 
global minimum of the effective potential appears at a positive value of a. 

There is freedom in the choice of the sign of the fermion mass term. In Fig. [2] we 
plot the behavior of the effective potential for rriu = —rn^. In this case, it depends on 
0"^ + (tt^)^ + (tt^)^ and (tt'^)^. In the figure, we plot the effective potential in the tt^ = 
plane and the a = tt^'^ = plane. It is seen that in this case, the effective potential exhibits 
behavior similar to that the sigma axis, vr" = 0, and the charged pion axis, cr = tt'^ = 0. As 
seen in Fig. [2l the global minimum of the effective potential lies in the a = tt^'^ = plane. 
This implies that the neutral pion field develops a non-vanishing expectation value, (tt^) ^ 0. 
It is found that pion condensation takes place for |mu + md| ^ 0.020 MeV. Here, we consider 
only the leading-order in the 1/N expansion. Higher order corrections can modify the small 
difference of the lines in Fig. [2j 

If |mu| + |md| is fixed at a realistic value (~ 9 MeV) and both m^^ and are positive, 
the global minimum of the effective potential lies in the tt" = plane. An enhancement 
of the symmetry breaking in the vr" = plane is mainly caused by the isospin symmetric 
terms, which are proportional to |mu + md\ in Eq. (12- 101) . Only a develops a non- vanishing 
expectation value. For small values of |mu-|-md|, the situation is different. The higher-order 
terms, specifically, O(m^), in Eq. (12- 101) break the isospin symmetry and yield non-trivial 
corrections to the effective potential. These corrections can move the global minimum of 
the effective potential away from the vr" = plane. If we take ^ —fna, the neutral 
pion develops a non-vanishing vacuum expectation value, and the vacuum state breaks the 
discrete symmetry under the parity transformation, as pointed out in Ref. 8). 

In such calculations, the mass scales of the effective potential are determined from the 
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-6.62x10-= - 

0.26 0.27 0.28 

a,7t"(GeV) 

Fig. 2. Behavior of the effective potential (j2-10p for mu + md = MeV in the vr^ = plane (dashed 
curve) and the a = vr^'^ = plane (solid curve). 



physical observables. For this purpose, we calculate the pion decay constant /tt and the pion 
mass rriTT and fix the three-momentum cutoff Af and the coupling constant G.^'^^ Here we 
set the current quark mass as + = 9 MeV and m = (mu + In this case, the 

minimum of the effective potential is at cr = (cr) and -k"" = 0. Taking the partial derivative 
of the generating functional (I2-6I) with respect to tt twice, we obtain the two-point vertex 
function for the pion field, r'i^\p). At a = (cr) and tt = 0, it is given by 



N 
G 



N 



Ac 



rtr 



G ^ 27r2 



2/lfv/^f2 + (a + m)2 



?2 >y 4((T + m)2 — p2 arctan 



^ylf2 + (a + m) V4(f^ + my - 



- {2{a + mf -p'^]\n 
for < < 4((T + mf 



ylf + y/ylf^ + (g + 
(T + m 



(2-14) 



The wave function renormalization, Z^r, is defined as ri'^\p) = ^p'^ + 0{p^). Evaluating 
the coefficient of in Eq. (12- 141) . we obtain 



Z-' 



N 
2^ 



In 



Ai+ v//lf2 + {a + 



m] 



a + m 



^A^ + (a + m)2 



The pion decay constant is given by 



(2-15) 



(2-16) 



We also use the following model-independent relationship, which is known as, the Gell- 
Mann-Oakes-Renner relation: 

A V = ^ma. (2-17) 
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Table II. Influence of the isospin breaking on the charged and neutral pion mass. 



(MeV) 


rrid 
(MeV) 


Min. oiVmia^n") 
(MeV) 


ma- 
(MeV) 


(MeV) 


(MeV) 


m^j-i — rrij^o 
(MeV) 


4.5 


4.5 


(o-) = 282.9 


590.4 


134.0 


134.0 





2.5 


6.5 


(a) ^ 282.9 


590.4 


134.0 


133.9 


0.1 


-4.5 


4.5 


(vr'') = 271.9 


8.948 


8.948 


538.9 


-530 


-2.5 


6.5 


(a) ^ 276.9 


562.4 


88.75 


88.29 


0.46 



Inserting the experimental values /,r = 91-9 MeV and tjIt^ = 135 MeV, which are neutral 
pion values, and m = 4.5 MeV into Eqs. (I2-16P and (12-171) and solving the gap equation, 
dVm/da = 0, we fix the cutoff scale and the coupling constant to A{ = 0.697 GeV and G = 
24.8 GeV~^. These parameter values were determined without QED corrections. However, 
the QED scale is much smaller than the typical scale for pions. It is also appropriate to use 
these parameter values in the gauged NJL model considered below. 

Using Eqs. (12-lOp and (I2-I6p . we can calculate the mass difference for the charged and 
neutral pions. 



2 - m^o^ = 2NZ^ 



9(7ri'2)2 5(vr3)2 



(2-18) 

(<7>,(7r<'> 



where we evaluate the right-hand side of this equation at the minimum of the effective 
potential, V^(cr, vr"). We list the pion mass and the mass difference in Table II. Due to the 
SU{2) isospin symmetry, the charged and neutral pions have the same mass for m„ = m^. 
In the case m„ = —2.5 MeV and = 6.5 MeV, the charged and neutral pion masses 
become smaller than in the positive up and down quark mass case. For = —4.5 MeV 
and = 4.5 MeV, the vacuum expectation value is in the direction. Non- vanishing 
expectation value is developed for (vr^). Hence the roles of a and vr^ are exchanged. 

Here we use the cutoff and the coupling constant for m„ + rrid = 9 MeV. As seen in Table 
UTt the pion masses are quite different from 135 MeV for mu = —2.5 and —4.5 MeV. In these 
cases, Eqs. (I2-15P and (I2-17P are not satisfied with = 135 MeV and f-„ = 91.9 MeV. In 
order to reproduce realistic values for m.,^ and f^^ in the case with negative m„, we should 
modify the condition (12-170 and use suitable values of Af and G. 

§3. Gauged NJL model 

Since the electric charges of the up and down quarks are different, an electromagnetic 
interaction breaks the SU{2) isospin symmetry. To elucidate the infiuence of this isospin 
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breaking, we introduce the QED correction. The QED interaction causes the derivative in 
the Lagrangian density (12- ip to be replaced with a covariant one. Thus, the Lagrangian 
density is extended to that of the gauged NJL model ^ 

A = ^ + leQA^) - M]i^ + [{i,^^^ + {^i^.r^i^f] , (3-1) 

where Q represents the charge of quark fields, Q = diag(2/3, —1/3). In this Lagrangian 

density, the isospin symmetry is broken, except along the direction. Since the QED 

corrections come from the internal photon lines at lowest order, we should also consider the 
free Lagrangian for photons, 

>Cph = —F.^F^'^ - ^{d,An\ (3-2) 

and evaluate the theory with the Lagrangian 

C = /:i + ^ph. (3-3) 

First, we evaluate the path integral for the quark fields. Employing the auxiliary field 
method applied in the previous section, we introduce the fields cr and vr and calculate the 
generating functional. Then, replacing the derivative in Eq. fl2-6p by a covariant one, we 
obtain 

Z = j VaVnVAexp i j d^xC^^. + /^^^ + 

-i In det [i^^{d^ + ieQA^) -M-a- iTsrV"] | 
The second line in Eq. f l3-4l) is expanded as 



(34) 



i In det [i'^^{d^ + ieQA^) -M-a- i-i^r^Ti"] 

oo 

= itr ln(i7^9/, -a- i-i^r'^-K") + tr ^ J„, (3-5) 

n=l 

where Jn is given by 

m \i'~i^^o^ — a — i'^^T°''K°- + le J 

Because the electric charge, e, and the current quark mass, M, are sufficiently small, we 
evaluate Eq. fl3-6p up to n = 2. At this order, the current quark mass and photon dependent 



*' As discussed in Ref. 12), this model has a large anomalous dimension, and the coupling constant is 
walking (not running). Hence, it may be possible to renormalize the theory even in four dimensions. 
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Fig. 3. The graphs expressing Eq. (j3-7p . The circles represent the quark propagators with the 
constituent quark mass, and the wavy hnes represent photons. 



parts are separate. The quark mass dependent part is identical to that obtained in the 
previous section. We calculate the photon dependent part using dimensional regularization 
and get 



trJ; 



ph 



d X 



d^p d^k 1 
L tr 

(27r)4 i(27r)^ -fpkP - a - i-f^r^n'' 



eofi'Q-fpA^'ip) 



d'^p 



I d'x 



A^ip)A,i-p) [{p^p''-g^V)nip') 



(3-7) 

Note that there is no term proportional to A{p)'^tt^7t^ in Eq. (13 -Zp . The vacuum self-energies 
for the photon field, UijP') and IlT^T^ijP'), are given by the diagrams in Fig. [3l Performing 
the momentum integration for the internal fermion lines, we get 

(3-8) 
(3-9) 



i7(/) 


Sao 
~ 277r 


'1 

e 


, 47ru2 8 

- In — hr + - 
ct'2 3 






_ 2ao 1 

TT ' 




f In + 2 





where the regularization parameter e is given by e = (4 — Z})/2, and we have defined 
eo/i" = e, ao = eo^/(47r) and /i = a/1 — 4(t'2/ (p^ -|- ley 

These self-energies are divergent in the four-dimensional limit. To obtain a finite result, 
we must renormalize the photon self-energy. Here we introduce the following counter-terms 
into Eqs. (lF7l) - (l3^ : 

1 



AC 



4 



(3-10) 



where the wave function renormalizations and Ztttt are defined by Z3 — 1 = —Nn{0) and 
= — (A^/2)i77r7r(0). We shift the pion fields as 7r^(a;) {t^^) +7r^(x)/-\/iV and normalize 
the photon field as N/?, ■ A^. Taking the large N limit, i.e. 1/N <^ cq, we obtain 

i j d^x{Cpi, + AC) - iNiijf - ^ y d^xAp {9'""d^ - d'^d") (l + Nn^{x)) 
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+-9^9" + Ng^-'TT+TT-nf^ix) 



A, 



where the renormahzed self-energies read 



i7^(p2) = n{p') - 77(0) 



ha 



3 2^ ' h-1 



and 



2a 



TT 



2-hln 



h + 1 
h-1 



(3-11) 

(3-12) 
(3-13) 



Performing the path integral over and the Wick rotation, t —2x4, we obtain 



j VAexp i j rfV(£ph + AC) - iNiijf 



exp 



iN 

— tr {ln(5^, + N-K+'K-nf^D^,) - \nD^,] 



(3-14) 



where tr represents the trace operation over the space-time coordinates and D^y represents 
the photon propagator 

1 



~k^ 



k^ky 

~k^' 



(3-15) 



Substituting Eqs. fl3-12p and (13-131) into Eq. (13-141) . we obtain the generating functional in 
the Landau gauge, 



Z ~ + / VcrV-K exp 



iN 



d X 



d^p 



X <3\n 



p' + 4(t'2 ( 1 + 



4aN\ AaN 



277T J 



3n 



-TT TT 



+ In 



p' + 4a" 1 + 



4aN\ 



21tx 



41n(p^ + 4a 



/2^ 



(3-16) 



In Eq. (13-161) . we omit terms which do not depend on a and tt and apply the following 
approximation: 



h(p^)\n^^^^l±^ 



1 + ^+ ' 



?)h{p'^Y 5h{p 



,2^4 



+ - - 



for p'^ ~ 0. 



After the momentum integration, the effective potential is found to be 

\/((T,7r°) = Kn(a,7r'^) + l^ga,(a,7r'^), 

where Vgs,^{a,7!'"-) is given by 



(3-17) 



(3-18) 



1 



192n' 



[3/ (4a'2(l + 4aN /27n) - (4aA^/37r)7r+7r"; V) 



+/ {4a'\l + 4aN/277c); V) - 4/(4^'^; A^' 



(3-19) 
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-6.2x10 




,-5 



> 

O -6.4x10" 



-5 



o 



> 



-6.6x10 



0.24 



0.26 



0.28 



0.3 



0.24 



0.26 



0.28 



0.3 



7t''(GeV) 



o(GeV) 



(a) 



(b) 



Fig. 4. Behavior of the effective potential (|3-18p for m-„ + md = 9 MeV. In (a), the black curve 
represents 1/(0, vr") at e = 0. The red and blue curves represent 1/(0, vr^'^, 0) and V{0,0,-k^), 
respectively. In (b), the solid and dotted curves represent y(cj, 0,0) at e = and y(cr, 0,0), 
respectively. 

Here, the function /(s^; t^) is defined in Eq. fl2-lip . and A{ and Ap are the three-momentum 
UV cutoffs for the quark and photon momenta, respectively. 

To study the QED corrections for the phase structure, we numerically evaluated the 
effective potential fl3-18p . In these calculations, we set A{ = Ap and chose the average of the 
current quark mass to be m = 4.5 MeV. The behavior of the effective potential is shown in 
Fig. m It is seen that the QED correction slightly suppresses the symmetry breaking for 
Via^TT" = 0) and V{a = 0,7t^^^ = 0, vr^) but enhances it for V{a = 0,7T^'^,n^ = 0). This 
result for a and ir^ is different from these obtained in previous studies. This result depends 
on the procedure used to regularize the theory. However, the global minimum of the effective 
potential lies on the line (vr") = 0. Only the scalar composite field a acquires a non- vanishing 
vacuum expectation value, and the orientation of the chiral symmetry breaking remains in 
the a direction even with the QED correction at the one loop level. Taking the massless 
limit, m — i> 0, the effective potential in the a direction coincides with that in the tt^ direction. 

Because of the QED correction, the second derivative of the effective potential in terms of 
TT^'^ becomes smaller than that of vr^. Thus, to leading order in the 1/N expansion, we obtain 
a result for the charged and the neutral pion mass difference whose sign is the opposite of 
the observed value. We consider this further in the next section. 



As is well-known, a pion-photon interaction plays an essential role in determining the 
mass difference between charged and neutral pions.^^-* Such an interaction is not included 



4. Charged and neutral pion mass 



13 



Fig. 5. Leading-order diagrams of the vacuum polarization for photons. The wavy hues and the 
dashed hne represent the photon and charged pion propagators, respectively. 



in the NJL model to leading order in the expansion. Here, we phenomenologically 
introduce higher-order corrections through the kinetic term for mesons,^''-* 



1 

2iV 



(4-1) 



where we fix = 3, because this term appears at the next-to-leading order in the 
expansion. The total Lagrangian is given by 



C — C{ + £ph + Cs- 
The expectation value of the charged pions induces the photon mass 



(4-2) 



e^Ti^T^ /3. This seems to suppress the QED correction for the effective potential 



V{a = 0,7r 



1,2 



0). 



To determine it, we calculate the effective potential for the theory with the Lagrangian 
(14-21) . The leading-order corrections to the vacuum polarization amplitude are represented 
by diagrams in Fig. O Only the charged pions contribute to these corrections. As shown 
in Refs. 15) and 16), the neutral pion mass should be invariant under the QED correction. 
After the normalization a/A^/3A^, we define the effective action by 



j DAexp i j {Cp^ + + AC) - iNtijf 



exp {iNr[a, tt"]} 



To leading order in the 1/A^ expansion, it reads 



r[cr,7r"l ~ -Indet 

L ' J g 



In the Landau gauge the effective potential is given by 

1 



Vgau(0-, Tt") 



2 / rf% L ' 



1 



12 J i(27r)' 



3 In 



4cr'^ I 1 + —r— I : TT+TT - 



+ ln 



4a'2 ( 1 + 



277r 



V 



Tin / 37r 



- 4 ln(4(T'^ - P ) + 3 ln(mA^ - V 



(4-3) 



(4-4) 



(4-5) 
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a,n (GeV) 

Fig. 6. Behavior of the effective potential ()4-6p for m-u + ?n.d 
V{a, vr" = 0), and the red and blue curves represent V{a ■ 
0, vr^), respectively. 



9 MeV. The black curve represents 



0) and V{a = 0,tt 



1,2 



Thus the total effective potential can be written 

V{a, 71-) = V^{a, + rgau+pi(t^, vr-^), (4-6) 
where Vgau+pi(o-, tt'') is given by 

V.+pi(^, vr-^) = ^gau(^, vr^) + ^/(^a'; Al). (4-7) 

We numerically analyzed the effective potential for Af = Ap at m = 4.5 MeV in this case 
also. The resulting behavior is plotted along the a, vr^'^ and vr^ axes in Fig. [61 Comparing 
these results with those in Fig. |H we find that the effective potential along the vr^'^ axis, 
V{a = 0, TT^'^, TT^ = 0), is higher. The curvature of the effective potential increases along the 
vr^'^ direction around the minimum. Thus the charged pion acquires the larger mass through 
the pion-photon interaction. Inserting Eq. (14 -Op into Eq. (12- 181) . the mass difference between 
the charged and neutral pions is found to be 

m^± - m^o ~ 7.9 MeV. (4-8) 

This value is larger than the observed one. We note that this result depends on the NJL pa- 
rameters Af and Ap. This implies a non-negligible uncertainty in our result. For this reason, 
if we consider phenomena in which the isospin breaking of the pion mass may contribute, 
we must include higher-order corrections in the 1/A^ expansion. 

§5. Gauged NJL model at finite temperature 

It is conjectured that the order parameters of the symmetry breaking (a) and (vr'*) disap- 
pear at the critical temperature in the chiral limit, M ^ 0. The broken chiral symmetry is 
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restored at higher temperature, and a state of deconfined partonic matter is reahzed. Here 
we consider a thermal system below the temperature at which the expectation values (a) 
and (vr") disappear. 

Following the standard procedure of the imaginary time formalism/^^'^^-' we introduce 
the temperature into our model. We then obtain the effective potential for the NJL model 
(12- II) to leading order in the expansion at finite temperature as^^-* 



(2^ 

2, . ^ r d^k 1 



71= — oo 

OO 



(mu + md)cr ^ / 



' ■■^"^^e.y (2vr)3 (a;F)2 + + ^,2 



71= — oo 



oo 



d^k 1 



X [{m^^ + m,') mr + k' + a'^} 

-2 {{mj + m^)a'' + (m^ - md)'7r+7r-}] , (5-1) 

where (5 is the inverse temperature, (5 = 1/T, and the discrete variable is given by 
— (^//^)(2'^+l), with anti-periodic boundary conditions for fermions. Using the formulae 



n=— oo 

oo 



In {(cij^)^ + x^} = 2 Incosh — X + [constant], (5-3) 



2 

n=— oo 

we carry out the summation in Eq. (15 -ip and get 



Vi{a,7r'^) = ^[a' + {nn']-^ / ' dk kHn cosh (^Vl^T^ 

-7r^{^n + md)a / dk - tanh -Vfc^ + 

27r^ Jq y/k^ + a'^ 2 

-(m,,^ + md^) / dk tanh — -[_ n-^S 

' io VFT^ 2 



47r' 

-Ai ^2 

X / (iA; 







A;2 + a'2 



1 . , /5 



VF+o^ 2 



tanh — Vft^+cr 



12 



- — sech ^ — v^P~+~a^ 
2 2 



(5-4) 



16 



o 

cog 

> -1x10-5 - 

-0.2 -0.1 0.1 0.2 

a(GeV) 

Fig. 7. Behavior of the effective potential (l5-4p . The red curve represents Vm{a,Tr'^ = 0) for 
m-u + TTT-d = MeV, the blue curve represents vr" = 0) for + = 4 MeV, and the 

black curve represents Vmic, vr" = 0) for m-u + = 9 MeV at T = 170 MeV. 



For mu = 2.5 MeV and = 6.5 MeV, the global minimum of the effective potential lies 
in the vr" = plane. In Fig. [71 we plot the behavior of the effective potential VJl{a, vr" = 0) in 
this plane. In the ground state, only the scalar field a develops a non-vanishing expectation 
value. As seen in Fig. [3, a small mass difference between up and down quarks induces a 
larger difference for the effective potential near the critical temperature. Since the effective 
potential in the a = plane, Vj^{(7 = 0, vr'^), has only a weak dependence on the current 
quark masses, and ma, Vj^{a = 0, vr*^) for + ma on the order of a few MeV has a 
shape that is similar to that for + = 0. In the limit m-u + rrid — > 0, isospin symmetry 
breaking appears in terms of O(m^). Thus the shape of the effective potential V^{a = 0, vr'*) 
is almost the same as that of Vj^{a, tt" = 0), which is plotted by the red curve in Fig. [3 The 
difference between the shapes of the effective potentials in the vr" = plane and the a = 
plane is extremely small. In fact, we cannot distinguish these shapes in Fig. [71 

It should be noted that the isospin symmetry breaking arises from the difference between 
the masses = 2.5 MeV and rrid = 6.5 MeV. This causes the small differences between 
(7r^'^)o.=o,7r3=o and (vr^)a.=o,7ri>2=o- The ratio for the local minimum. 



(5-5) 



is approximately 0.05% for = 2.5 MeV and = 6.5 MeV at T = 170 MeV. The 
effective potential at the point cr = 0, vr^'^ = 0, vr^ = (7r^)o-=o,7ri.2=o is slightly smaller than 
that at o" = 0, vr^'^ = (7r^'^)o-=o,7r3=0) ^'^ = 0. The ratio for the depth of the effective potential, 
rv given in ( 12-13p . is approximately 0.2% at T = 170 MeV. In the case + rrid < 0.0038 
MeV, we again observe neutral pion condensation. 

In the gauged NJL model, the situation is somewhat complicated. Due to the Debye 
screening photon mass mp'^ = (5/54) e^T^ + 0(T), the photon propagator is modified at 
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finite temperature. In the case of electrostatic shielding, po — > 0, it is expressed as 

n - pT I ^ I ^ ^^^^ rfi.fil 

+ Inip'^ ^ p^ p^ 

Pl^Q, Pl^^O, (5-7) 

- V - ^ - (5-8) 

Replacing the first two terms of Eq. f l4-4p with Eq. (15-60 . we obtain the QED corrections 
to the effective potential at finite temperature in the Landau gauge, 

n= — OO " \ ' L V y 

-3 In + p2 + ^ 2 In + + mA^} 

+ ln{{u;^f + p^ + mA^ + 2mp^}], (5-9) 

where the square of the four momentum p"^ is written p"^ = (uj^)"^ + P^^ discrete 
variable is given by uj^ = {tt/ (3)2n, from the periodic boundary conditions for bosons. 
From the above, the effective potential of the gauged NJL model is found to be 

V^<y, tt") = Vl{a, vr") + vr'^). (5-10) 



The photon dependent part Vlau+pi(o", vr") is given by 



gau+pi 



i 



3 insmh — a / + 4(t'^ tt+tt" 

2 V Stt 



+ In sinh — ^Jp^ + m\ + 2m/3 



—3 In sinh — \/ p^ + 4cr'^ + 2 In sinh — y + 



(5-11) 



where we have used the formula 



°° 3 

^ In {(cj^)^ + x^} = 2 In sinh + [constant]. (5-12) 

n=— OO 

We numerically calculated Eq. (15 -101) with Eqs. fl5-4p and (15- 111) just below the critical 
temperature. In Fig. [HI we restrict the parameter space to the region in which two of the 
composite fields cr, tt^'^ and vr^ vanish and plot the effective potential fl5-10p for T = 170 
MeV. For + md = 9 MeV, we find the global minimum of the effective potential at 
= 0. Only the composite field a condenses even at finite temperature. The 
orientation of the chiral symmetry breaking is again preserved under the QED correction. 
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a,n(GeW) 

Fig. 8. Behavior of the effective potential ([500]) for m^ + md = 9 MeV at T = 170 MeV. The black 
curve represents V^{a,TT"' = 0), and the red and blue curves represent V^{a = 0, vr^'^jVr^ = 0) 
and V^{a = 0, tt^'^ = 0, tt^), respectively. 
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Fig. 9. Behavior of the global and local minima of the effective potential (j5-10p for + = 9 
MeV. The black curve represents the location of the global minimum, ((T),r°=o- The red and 
blue curves represent the local minima (7r^'^)^=.^3=o ^^'^ {'^^)a=n'^-'^=Oj respectively. 

In Fig.[9l we plot the minimum of the effective potential in the restricted parameter space, 
V^{a,n'' = 0), V^{a = 0,TT^'^,n^ = 0) and V^{(r = 0, tt^'^ = Cvr^). As is clearly seen in 
the figure, the expectation value (cr) disappears smoothly at higher temperature. Thus, the 
transition from the broken phase to the symmetric phase is crossover. This is a well-known 
feature of the theory with a current quark mass. It should be noted that the pion fields have 
vanishing expectation values in the ground state for m^^ + m^ = 9MeV. It is conjectured that 
the chiral symmetry is restored in the early universe. Non-vanishing expectation values of 
the pion fields (7r^'^)o-=o,7r3=o and (7r'^)o-=o,7ri.2=o ^lay be realized during the transition period 
of the chiral symmetry breaking, i.e. "disoriented chiral condensation" .^^^'^^^ 

Near the critical temperature we observe that the ratio ry, which is defined in Eq. f l2-13p . 
is much larger than in the T = case. It is believed that the influence of the isospin breaking 
may be observed in phenomena near the critical temperature. 

Inserting Eq. (15-101) into Eq. (12-181) and using Eq. (12-151) . we evaluate the mass for the 
TT meson or a soft mode which corresponds to tt. In Fig. [10] we plot the behavior of the 
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(a) Charged and neutral pions masses. (b) Pion mass difference. 

Fig. 10. Behavior of the pion mass for rriu + ma = 9 MeV. The red and blue curves represent the 
charged and neutral pion masses, respectively 

pion and sigma masses and the pion mass difference as functions of temperature. As seen in 
Fig. 10(a), the pion and sigma masses exhibit behavior similar to that found in the previous 
studies Refs. 3) and 23). A minimum of the sigma mass is found near the temperature 
T = 200 MeV. We call this the critical temperature Tc- Above Tc, a pion can exist as a soft 
mode.^'*-' The QED correction enhances the charged pion mass, even at finite temperature. 
The pion mass difference has a maximum value below the critical temperature. 

§6. Pion condensation and dynamical CP violation 

In the present section, we consider the more precise structure of the pion condensation 
and discuss the dynamical CP violation induced in the four-fermion interaction model. 

The QCD Lagrangian is invariant under the CP transformation except for quark mixing 
terms. Here we consider the NJL model (12- ip as a low energy effective theory of QCD. 
This model does not contain any quark mixing terms. Thus the Lagrangian density (12- ip 
is CP-invariant. Because a pseudo-scalar operator is CP odd, the expectation value of the 
pion field transforms under CP as 

(tt") ^ - (|^^^75r>^ ^ - (tt^) ^ ^^V^nsr^^ . (6-1) 

Thus, the non-vanishing expectation value of the pion field spontaneously violates the CP 
invariance. Such a mechanism for violating the CP invariance was suggested long ago by 
Dashen.^^ A simple four-fermion model to realize such a mechanism was proposed in other 
contexts in Refs. 25) and 26). 

As shown in the previous sections, pion condensation is observed for a vanishing average 
of the current quark mass, + = 0. In Table III we list the values of the effective 
potential at the minima on the a, tt^'^ and tt^ axes for + = at T = and T = 170 
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Fig. 11. Phase structure in the m-T plane. 

MeV. The global minimum lies on the (vr^) axis. The temperature enhances the ratio 
defined by 

_ V{{a) = 0, (7r"))-V-((a),(7r-) = 0) 

V{{a) = 0,{n-)) + V{{a),{n'^) = Oy ^ ' 

Hence we see that a non-trivial CP violating phase is realized for the NJL model. 

We also evaluated the effective potential for small but finite values of |mu + md| (= 2m) 
and obtained the phase structure of the chiral symmetry breaking in the m — T plane. As 
seen in Fig. [Ill the CP violating phase, (vr^) ^ 0, appears within a narrow parameter range. 
In Fig. [TT] we cannot observe the restoration of the broken chiral symmetry, except for 
mu + rrid = 0. The expectation value (cr) disappears smoothly at higher temperature, e.g. 
in the case depicted in Fig. [91 

Because our original Lagrangian has no ^-term, the dynamical CP violation takes place 
in the NJL model with two fiavors of light quarks. The CP violating phase is observed for 
|mu + m^l < 0.020 MeV at T = 0. Below, we consider two assumptions concerning the 
current quark mass. 

With the assumption |mu + md| < 0.020 MeV and \m^\ + |md| > 0.020 MeV, the up and 
down quark masses have opposite sign. The sign of the quark mass is reversed by the discrete 
chiral transformation u —>■ 75M. Thus, the theory with a mass term, |mu + md| ^ 0.020 MeV, 
is obtained through the chiral transformation. Under this transformation the expectation 
value of the pion field vanishes, and therefore the ground state preserves the CP invariance. 



Table III. The depth of the effective potential and the ratio for m-u + ?Tid = 0. 



T 


V{{a),0,0) 


m (7ri'2),0) 


F(0,0,(O) 


^1,2 


ra 


(MeV) 


(GeV^) 


(GeV) 


(GeV^) 









-6.60 X 10-5 


-6.60 X 10-5 


-6.61 X 10-5 





4.24 X lO-'^ 


170 


-1.30 X 10-'^ 


-1.30 X 10-^ 


-1.33 X 10-^ 





0.00943 
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However, the discrete chiral transformation is anomalous. It produces the CP violating 
6'-term with the angle tt in QCD. This is known as a dynamical mechanism violating CP 
invar iance.^-* This term is inconsistent with the experimental upper bound for the neutron 
electric dipole moment. We thus conclude that another contribution to the neutron electric 
dipole moment, for example contributions from heavier particles, should be introduced. 

With the other assumption, |m„ + m^il < 0.020 MeV and < \m^\ + |md| < 0.020 MeV, 
we cannot obtain the Lagrangian with |mu + m^l > 0.020 MeV through any redefinition of 
the quark fields. In particular, if either the up or down quark is massless, any ^^-term can 
be rotated away. This case has been investigated in Refs. 27) and 8). With the assumption 
|mu| + |md| < 0.020 MeV, the model is entirely incapable of accounting for the real meson 
properties. But, nevertheless, it is interesting to suggest the possibility that dynamical CP 
violation occurs at low temperature without introducing the strong CP problem. 

§7. Conclusion 

We have investigated the symmetry properties of the NJL model with two flavors of 
quarks. The exphcit breaking of the global flavor symmetry SUi,{2)® 5'C/r(2) was introduced 
from the current quark mass and the QED interaction. We have considered the explicit 
breaking of this symmetry up to O(m^) and O(e^) and calculated the effective potential to 
leading order in the 1/N expansion. 

Evaluating the effective potential, we have determined the precise structure of the chiral 
symmetry breaking. The effective potential has terms proportional to the average of the 
current quark mass at 0{m). We found that a positive quark mass enhances the chiral 
symmetry breaking for a positive a and suppresses it for a negative a. A negative quark 
mass has the opposite effect. For a realistic quark mass, mu+md 9 MeV, it was found that 
the ground state is characterized by tt" = 0. Only the scalar composite state a develops a 
non-vanishing expectation value. This situation does not change when we introduce the one- 
loop QED corrections and the finite temperature corrections. Therefore we conclude that 
the orientation of the chiral symmetry breaking is invariant for a sufficiently large magnitude 
of the average of the quark mass, |mu -|- md| > 0.020 MeV. 

The difference between the current quark masses of the up and down quarks begins to 
affect the effective potential at second order, 0(|mu— mdp). It contributes to the tt^ direction 
and causes isospin breaking for the effective potential, but the ratio of V{{tx^''^)) to V^((7r"^)) 
is extremely small, 0.008% at T = and 0.2% for = 2.5 MeV and = 6.5 MeV at 
T = 170 MeV. 

These effects cannot change the ground state, but there is a possibility that the state 
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passes through a local minimum of the effective potential at the time of the symmetry 
breaking in the early universe and the pion fields thus temporarily have a non-vanishing 
expectation value. Highly excited states arc produced near the critical temperature; these 
excited states may move from the ground state with positive {a) to a finite (tt) state. This 
phenomenon is called "disoriented chiral condensation". ^^''^^^ 

To explain the charged and neutral pion mass difference, we have to consider higher-order 
contributions in the 1/N expansion. Here we phenomenologically introduced a pion kinetic 
term. Then we obtained the correct sign for m^± — m^o . The value of m^± — m^o increases 
near the critical temperature. The neutral pion direction is more stable than the charged 
pion direction, even at finite temperature. It is conjectured that the production number of 
neutral pions is larger than that of charged pions when the QCD ground state passes from 
the quark ghion state to the hadron state. 

If we assume that the average of the quark mass is sufficiently small, the enhancement 
of the chiral symmetry breaking in the a direction will not be large. We cannot ignore the 
isospin breaking terms that come from the mass difference even in the ground state. Neutral 
pion condensation can be realized in a restricted parameter range. We can consider the 
case in which only a develops a nou- vanishing expectation vahic through the chiral SU{2) 
transformation. Then, the only difference is found in the ^^-term. Thus the pion condensation 
studied here affects only phenomena in which the ^7^(1) anomaly contributes. Such pion 
condensation may not be realized in the present universe, but it would be interesting to 
study the dynamical origin of CP violation and the critical phenomena in the early universe. 
Especially for |mu| + \md\ < 0.020 MeV, we find the possibility of dynamical CP violation. 

There is a possibility that the composite operator tt'^tt" also develops a non-vanishing 
expectation value. In the leading-order analysis of the expansion we cannot evaluate 
the condensation (tt^tt^). Thus, an extension of the model or the analysis is necessary to 
elucidate such a phenomenon caused by the pion dynamics. 
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